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Abstract. A class of left-invariant second order reversible sys¬ 
tems with functional parameter is introduced which exhibits the 
phenomenon of robust integrability: an open and dense subset of 
the phase space is filled with invariant tori carrying quasi-periodic 
motions, and this behavior persists under perturbations within the 
class. 

Real-analytic volume preserving systems are found in this class 
which have positive Lyapunov exponents on an open subset, and 
the complement filled with invariant tori. 


1. Introduction 

We study a family of second order dynamical systems on a locally 
homogeneous Riemannian space M, modeled on a special solvable Lie 
group. The simplest example is the geodesic flow of a left-invariant 
metric. Our class generalizes the examples discovered by Butler [ Bui] , 
and Bolsinov and Taimanov IEQ. In these examples the complete 
integrability of the geodesic flow in the tangent bundle TM is accom¬ 
panied by highly non-integrable behavior on an invariant submanifold 
of codimension n = dim M. The dynamics there is the suspension 
of a toral automorphism, and in |B-Tj the hyperbolic automorphism 
is chosen, which leads to an Anosov flow. The presence of an Anosov 
flow as a subsystem guarantees the positivity of topological entropy. 

We show that such a behavior extends to a larger class of left- 
invariant second order systems. This class is parametrized by a matrix 
L and a functional parameter F: a smooth vector field in the unit ball 
of the Euclidean space. We call them L — F systems. The crucial prop¬ 
erty of L — F systems is their J-reversibility, where J is an appropriate 
involution of the tangent bundle TM. Let us recall that a system is J- 
reversible if the involution J conjugates the forward in time dynamics 
with the backward in time dynamics. 

There is a vast literature devoted to J-reversible systems. The survey 
paper of Lamb and Roberts ra contains an extensive bibliography. 

In particular there is a version of the KAM theory for the J-reversible 
systems. It goes back to Moser [Mj . and Sevryuk [5]. In our case we 
establish robust integrability: it persists under any small perturbation 
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as long as we stay in the family of L — F systems. Note that this family 
is parametrized by an infinite dimensional Banach space of vector fields 
F. What is notable is that we do not assume volume preservation, 
the symmetries imposed on the system force the integrability, and the 
occurrence of a finite absolutely continuous invariant measure. This 
measure has a density with respect to the Liouville volume which is 
only C°°, and typically no real-analytic invariant density exists (Section 
8 ). 

The ,/-reversible KAM theory would give us large subsets of quasi- 
periodic motions for perturbations which are not left-invariant, as long 
as they are /-reversible. We were unable to check the non-degeneracy 
of the unperturbed system required for the application of the KAM 
theory. However we conjecture that the non-degeneracy does hold for 
most systems under consideration. 

Butler, |Bu2]. |Bu3]. used the mechanism discovered in EZEI to ob¬ 
tain C°° examples of integrable volume preserving systems with posi¬ 
tive metric entropy. In our class we find whole families of real-analytic 
systems with positive metric entropy and a subset filled with quasi- 
periodic motions, open but not dense (Section 9). 

The phenomenon of robust integrability, accompanied by positive 
topological entropy occurs already for geodesic flows of linear connec¬ 
tions. The generalization of the geodesic flow of the Levi-Civita connec¬ 
tion to more general linear connections was discussed in |P-Wj . Such a 
generalization appears naturally in the study of Gaussian thermostats, 
a class of systems introduced by Hoover [H]. The paper of Gallavotti 
and Ruclle [G-R] introduces the Gaussian thermostats in the physi¬ 
cal context. In particular in our class of systems we find a Gaussian 
thermostat with the following paradoxical behavior (Section 10). For 
small kinetic energy the system is asymptotic to an Anosov flow, which 
has a large codimension in the whole phase space. For larger values of 
the kinetic energy the system undergoes a drastic change, it becomes 
integrable: an open and dense subset in the phase space is filled with 
quasi-periodic motions. The Anosov subsystem is still present, but for 
a subset of initial conditions of full Lebesgue measure the solutions stay 
away from that chaotic subsystem, and fill densely invariant tori. 

The author would like to thank Gabricla Tereszkiewicz, Wojciech 
Czernous and Adam Doliwa for very useful discussions. Preliminary 
results of this work were presented at the seminar of Lai-Sang Young 
at the Courant Institute in November 2014. We are grateful for her 
insights and hospitality. 

2. The configuration space 

Our configuration space is a locally homogeneous Riemannian space 
modeled on a special Lie group G. We start its description with the 
Lie algebra g. We assume that dim g = n+ 1 and that g contains an n 
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dimensional abelian ideal go- We choose an arbitrary scalar product in 
g, and let b denote a unit vector orthogonal to g 0 . Since the ideal g 0 is 
assumed to be abelian the Jacobi identity imposes no conditions on the 
operator L : g 0 —> g 0 , L = ad &. At this stage we place no restrictions 
on the operator L. Later on we will consider various special cases. We 
endow the Lie group with the left invariant metric determined by our 
choice of the scalar product in g. 

The Levi-Civita connection V on the Riemannian manifold G can be 
expressed as a tensor on g. It can be calculated directly, which is done 
in the fundamental paper of Milnor jMij , where extensive explanations 
can be found. The formulas read 

Vb& = 0, V 6 f = for f eg 0 , 

' 11 V ? 6 = , V^C = <S£, C )b, for £, C G g 0 

where S = + L*) and A — |(L — L*) denote the symmetric and 

skew-symmetric parts of the operator L. 

The Lie algebra has an important automorphism K : g —>• g, where 
—A' is equal to the euclidean reflection in g 0 . There are only few Lie 
algebras with this kind of additional symmetry. It is not difficult to 
enumerate all of them. Our class of Lie algebras is singled out by the 
additional property that the invariant subspace of K is a subalgebra. 

Since the automorphism K is orthogonal it generates the automor¬ 
phism % of the Lie group G which is an isometry. This isometry will 
play crucial role in our discussion. Let us note that both K and % are 
involutive, i.e., K = K~ 1 ,% = X -1 . 

The Lie group has the following matrix representation, which we will 
also denote by G. It consists of matrices with the block form 

1 U L , w G R n ,w G R. 


The abelian ideal go corresponds to the normal abelian subgroup Go 
consisting of the matrices with u — 0. We obtain convenient coordi¬ 
nates (w, u ) G M" x R in the group G, which is the semi-direct product 
Go x Gi, of the abelian additive groups Go = R n and G\ = R. 

The Lie algebra of our matrix group consists of matrices of the form 


0 0 
f VL 


£ e R n ,7] e R. 


We will be using the linear coordinates (£, rf) G R n x R in the Lie 
algebra g. 

Let us choose a lattice To in Go of rank n, so that To\Go is an n- 
dimensional torus. Our basic configuration space M is the Riemannian 
manifold M = To\G. The manifold M = To\ (Go x G\) is canonically 
diffeomorphic to (To\Go) x G\ = T n x R. In particular each torus in 
this product has a canonical affine structure inherited from the abelian 
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subgroup Go- This foliation into tori will play an important role. We 
will refer to the leaves of this foliation as toral leaves. 

In general the Riemannian manifold M is only locally homogeneous, 
since the left translations do not factor to the coset space T 0 \G. How¬ 
ever the left translations by elements from the subgroup Go do factor 
onto M, and these isometries preserve all the toral leaves, acting on 
them as translations. 

The configuration space M = To\G is non-compact. In some cases 
the discrete subgroup T 0 can be enlarged to a discrete subgroup Y so 
that the resulting quotient space N = Y\M is compact. It is so when 
the one parameter subgroup e uL ,u G M, contains an automorphism of 
T 0 C R n . In general the necessary condition for the existence of a 
quotient with finite volume is the unimodularity of the group G, [Mi j. 
In our case the group G is unimodular if and only if tr L — 0. 

Let us assume for simplicity that A = e L is an automorphism of the 
discrete subgroup T 0 (and the torus T 0 \Go). The discrete subgroup T 
is generated by T 0 and A G G {A — e L can be considered as an element 
of the subgroup Gp w — 0,u = 1.) We will consider such subgroups 
and the resulting compact locally homogeneous spaces N = h\G. 

3. Second order left invariant equations on Lie groups 

We will need some general properties of second order equations on 
Lie groups, which are invariant under left translations. 

A second order equation on a manifold M is a dynamical system, 
i.e., a continuous flow x Y t ,t G R, on the tangent bundle TM, such that 
for a trajectory (x(t),v(t)) = TTx(O),u(0)), t G M, where x(t) G M 
and v(t) G T x(t) M, we have § = v(t). 

Let further G be a Lie group. The tangent bundle TG of the group 
is canonically diffeomorphic to the cartesian product G x g by the use 
of left translations. With this identification the derivative of a left 
translation on G, as a mapping of TG = G x g, is equal to the identity 
in the second factor. 

A second order equation on a Lie group is left-invariant if left trans¬ 
lations on the tangent bundle TG = G x g commute with the flow. 
In other words the trajectories of the flow are taken to trajectories 
by left translations: for any g G G and any trajectory (x(t),v(t)) = 
v L < (a;(0),u(0)), t G M, the parametrized curve (gx(t),v(t)) in G x g is 
also a trajectory of the flow. 

The structure of a left-invariant second order equation is described 
in the following 

Proposition 3.1. If a continuous flow T 4 on the tangent bundle TG = 
Gxg of a Lie group is a left-invariant second order equation then it is a 
group extension of a flow on g, i.e., there is a continuous flow ^ : g — * 
g and a continuous cocycle over ip* with values in G, h : M x g —> G, 
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such that 

i & t (9,v) = (gh(t,v),‘il> t (v)), for (, g,v)eGxg. 

Proof. Let ttg ■ TG —> G and 7r g : TG —* g be the projections associ¬ 
ated with the identification of TG with Gxg. We define the flow 
in the Lie algebra by 

il> t (v) = 'ir B (V t (g,v)) for any g e G. 

The mapping if* : g — » g is well defined because of left invariance, and 
hence it must be a flow. 

Further let 

h(t,v) = g-'irc (^(giv)) for any g e G. 

Again since it is well defined the cocyle property can be easily checked. 

□ 

The flow on g described in Proposition 3.1 will be referred to as the 
Euler flow. This choice of terminology comes from the Euler equation 
of the rigid body dynamics, see [A], Appendix 2. 

It follows from this Proposition 3.1 that if the Euler flow has a peri¬ 
odic trajectory through v 0 G g of period T then the map preserves 
the set G x {u 0 } and it is equal there to a right translation on G. 

4. The dynamical systems: the L — F systems 

The simplest dynamical system that is of interest to us is the geodesic 
flow, which we consider in the tangent bundle TM, rather than the 
cotangent bundle. Using the Riemannian metric we can identify the 
tangent and cotangent bundles, and the tangent bundle acquires the 
natural symplectic form cu. 

The isometry X : G —* G projects to M — To\G, and we denote it 
again as X. This isometry is involutive, i.e., X = AT 1 . 

The derivative DX :TM —>• TM is also involutive and it commutes 
with the geodesic flow : TM -> TM , tel, 

DX o T 4 = T* o DX, t e M. 

Let J \TM —tTM be the involution which is identity in the base M 
(it does not move points in M) and it is equal to minus identity in every 
tangent space. It is well known that the geodesic flow is J-reversible, 
i.e., 

J O T* = T _t O J, i el 

The involutions J and DX commute, hence their composition is also 
an involution, and we denote it by J = J o DX = DX o J. Since the 
geodesic flow commutes with DX it must be also J-reversible. 

The J-reversibility is the fundamental self symmetry that is present 
in our family of dynamical systems on SM. 
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The first generalization is from the Levi-Civita connection of a left 
invariant metric on G to a left invariant linear connection V on G. 
Such a connection differs from the Levi-Civita connection by a tensor 
B in g 

V a -F - VaF = B(X, F), X, F G 0. 

We assume that the connection V has two additional properties. Firstly 
we require that the parametrization of its geodesics is proportional 
to the arc length. We will call such connections para-metric , since 
they generalize the concept of metric connections , which have iso¬ 
metric parallel transport. A connection is para-metric if and only if 
(B(X,X),X) = 0 for every X G g. The geodesic flow of a para-metric 
connection V preserves the unit sphere bundle, and it will be denoted 
again by : SM — > SM. A discussion of geodesic flows of linear 
connections can be found in jP-W] . 

Secondly we assume that the symmetric part of the tensor B = V—V 
is invariant under the isometric involution K : g —> g, i.e., for every 
XGg 

B(KX, KX) = KB(X, X). 

If such a property holds we say that the connection is weakly DC- 
invariant. It follows that for a left-invariant and weakly DC-invariant 
connection the geodesic flow T* is again J-reversible. 

The equations of the geodesic flow can be written as 


dx 

dt 


v , V v v = 0, 


where x{t) G G is a parametrized geodesic. 

With the identification of the tangent bundle TG with G x g by left 
translations we get there the coordinates (w,u\£,r}) G M 2n+2 , which 
were introduced in Section 1. In these coordinates the above involutions 
are given by 


DX(w , u] f, rj) = ( -w , u; 7]), J(w, u; f, rj) = (w, u; -rj), 


We need to establish the form of a left-invariant, para-metric, weakly 
DC-invariant connection on G. Recall that 6 G 0 is a unit vector 
orthogonal to go- An arbitrary element A" G g can be written as 
X = £ + qb, £ G g 0 - 

Proposition 4.1. A connection V = V + B on the group G is left- 
invariant, para-metric and weakly X-invariant if and only if there is a 
linear operator C : g 0 — > 0o such that for any X — £ + qb, £ G 0o. 


B(X,X)=r i CZ-(CU)b 
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Proof. Without loss of generality we can assume that B is symmetric. 
We have for every X = £ + rjb, £ € go, 

0 = (B(X, X),X) = <BR, 0,0 + (), b) 

+ + 2;f(B((,b),b) + ,, 2 (B(b.b),0 + V 3 {B(b,b),b). 

Since this cubic polynomial in g must vanish we get for every £ £ go 
<£(£, £),&> + 2<£(£, b), £> = o, 2(£?(£, b), b) + (B(b, &),£) = 0, 
(£(£,£),£> = 0, (B(6,6), 6) = 0. 

Since the connection is weakly X-invariant we must have KB(b,b ) = 
B(b,b) and AT?(£,£) = B (£,£), and so we obtain further that 
B(b,b ) = 0 and .£?(£,£) _L g 0 . If follows readily that 

<£(£,&),&> = 0 , 

£(£,£) = - 2 <£(£, &),£>&. 

Putting C£ = 2£>(£,&) we get the desired formula. □ 

Using Proposition 4.1 and the formulas (1) we obtain by direct cal¬ 
culations the equations of the geodesic flow for our special connections. 


Theorem 4.2. For any para-metric left-invariant and weakly %-invariant 
connection, the geodesic equations in the coordinates (w, u; £, rj) in the 
tangent bundle TG are 


( 2 ) 


df 

dt 

dw 

dt 


VF(0, 
- e ul f, 


drj 

dt 

du 

dt 


V- 


where F(£) = L*f — Cf, and the matrix C depends only on the tensor 
B = V — V ; namely for X = £ + 776 , £ e g 0 we have B(X,X) = 
vCf-(Cf,0b.D 


Note that in the last Theorem any matrix C can occur with an 
appropriate choice of the connection V. 

The equations (2) factor to the Lie algebra g as the Euler equations , 
the first line of (2). Let us recall that the Euler equations are obtained 
by left translations of velocities along geodesics. It follows from the left 
invariance of the connection that the resulting curves in g must satisfy 
the Euler equation, see [X] , Appendix 2. 

Our final generalization is to replace the linear vector held F(£) = 
E£,£ £ g 0 = in ( 2 ) by a general (non-linear) vector held 

F ■ go —;> go- Such a dynamical system will be called an L — F 
system. In the special case of a left invariant para-metric connection, 
when the vector held F is linear, we will call it a quadratic L — F 
system. 

To summarize: the operator L : g 0 —* g 0 in the euclidean space go 
determines the Lie group G with a chosen left invariant metric. The 
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vector field F determines then the second order equations (2) in TG 
which are preserved by any left translation of G, and hence project 
naturally to TM, or SM. Note that to define an L — F system on SM 
it is enough to have the vector field F defined in the closed unit ball 
of 0 o- 

The crucial property of general L — F systems on SM (or TM) is 
their J-reversibility. 

Actually the L — F systems can be characterized as smooth second 
order left invariant equations on SM which are J-reversible. 

Proposition 4.3. Any second order left-invariant equations on SM 
which are also J-reversible define an L — F system. 

Proof. By Proposition 3.1 the equations factor from SM to the unit 
sphere §> n = {(£, 77 ) G g |£ 2 + rfl = 1}. The resulting Euler system is 
also J-reversible, where J(£, 77 ) = (£, — 77 ). 

Using the coordinates (£ 2 , • • •, £n, v) 011 the unit sphere in the neigh¬ 
borhood of the equator {77 = 0} we consider the smooth function 

u ^ = ff 

By the J-reversibility we obtain that the function U is odd in the 77 
variable, i.e., [/(£, — 77 ) = — TJ(£, 77 ). It follows that if U is a smooth 
function on the unit sphere then the function 

m = -7=' m, %/u r P) 

is well defined and smooth in the closed unit ball in g 0 , and £/(£, 77 ) = 

ri F(0- □ 


5. Periodic solutions in J-reversible systems 

For J-reversible systems there is a very convenient way of searching 
for periodic solutions. It goes so far back that it is by now a part of the 
mathematical folklore. It was formulated explicitly by DeVogelaere 
[HeVj . and Devaney El- We do it again for the convenience of the 
reader. 

Let us consider the subset T of fixed points of the involution J, 
{ p \ J(p) = p }. 

Theorem 5.1. Any trajectory of a J-reversible flow T 4 which visits T 
twice must be periodic. Moreover such a trajectory is invariant under 
J with the reversal of time. 

Proof. Let p 0 G T be such that there is t 0 > 0 with G T. We 

have 


T to (p 0 ) = J (^°(p 0 )) = V ~ t0 J(po) = *- to {Po). 
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Hence the trajectory p(t) = 'J /< (po) is periodic with the period T = 2f 0 . 
Moreover 

Jp{t) = J¥(p 0 ) = V-Vpo) = ^~\po) = 

□ 

The minimal period of a trajectory of po in the above proof is T = 2f 0 
if and only if f 0 > 0 is the time of the first return of p 0 to the set T of 
fixed points of J. 

We apply this principle not to the full L — F system but only to its 
factor, the Euler equations 

( 3 ) | = = 

The involution J descends naturally to the phase space of (3) (£, 77) G g, 
and we denote it again by J, J(f,p) = (£,— 77 ). Clearly the Euler 
equation is J-reversible. The set of fixed points of J is equal to T = 
0 o = {{^,V)\V = 01 - 

Let us consider the open unit ball B C go, with the boundary, the 
unit sphere, S = dB. 

Definition 5.1. For a smooth vector held F = F(() defined on the 
closed unit ball in g 0 , we say that a point Co hi the open unit ball B is 
escaping if the integral curve ( = ((s) of F through Co is defined in a 
finite closed interval [s_,s+] 3 0, C(0) = Co, C( s ) G B for s G (s_,s+), 
the endpoints C( s ±) belong to the unit sphere S = dB, and the vector 
held F is transversal to the unit sphere S at the endpoints of the 
integral curve, i.e., 

( 4 ) <nc(s-)), C(O) < 0, <F(c( a+ )), C(*+)> > 0. 

The integral curve through an escaping point is called an escaping 
trajectory. 

Let us note that for a given smooth vector held F the set of escaping 
points in the unit ball is open. 

We have the following crucial 

Theorem 5.2. For any escaping trajectory ((s),s G [s_,s + ] of the 
vector field f = F( C), C £ 0 o> there are periodic functions 17 (f) and u(t) 
with the period T = 21 0 such that 77 ( 0 ) = 77 (f 0 ) = 0, tt(0) = s_,u(f 0 ) = 
s + ,s_ < 77(f) < s + , and (C(f), 77 (f)) is a T-periodic solution of the Euler 
equation (3), where £(f) = C (77(f)). Moreover pit) is an odd function 
77( —f) = — 77 (f), f G M, and u(t ) is an even function 

u(t ) = s_ + / 77 (s)ds. 

Jo 

Proof. We restrict the Euler equation (3) to the upper half of the unit 
sphere in g, i.e., to the subset {(C,T 7 )|C 2 + V 2 — l/h > 0}. In this 
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submanifold we can use ( as coordinates, and we introduce there the 
time change ^ = 77 . After this time change the Euler equation becomes 
the following system 

( 5 ) t = no ■ 

as 

Hence the trajectory £(s),s_ < u < s + , of this vector field in the 
open ball B C 0 o gives rise to the trajectory (£(f), 77(f)), f_ < t < t + , 
of the Euler equation in the upper half of the unit sphere in g, with 
linp_> t _ 77(f) = lim,^ t+ 77(f) = 0. It may happen (and it does) that 
f± = ±00. This is excluded in our case by the condition ( 4 ) which gives 
us linp_^_ ^ > 0, lim t _^ + < 0. Hence we have a finite time interval 
[f_,f + ] and we can shift it to the time interval [ 0 , fo], fo = f+ — f-. 

By Theorem 5.1 this trajectory extends to the T = 2fo periodic 
solution of the Euler equation (^(f), 77 (f)) ,f £ R. Putting 77(f) = s_ + 
f*ij(s)ds we have £(f) = £( u(t )) for 0 < f < t 0 . Since the periodic 
solution is invariant under J with the reversal of time, we obtain that 
£(—f) = £(t), 77 (—f) = — 77 (f). It follows that J i f > tQ ‘ri(s)ds = 0 , and 
consequently 77(f) is a T periodic function 0 < 77(f) < Uq. Hence we 
have also £(f) = £( u(t )) for any f eR. □ 

In simple terms what this proof reveals is that an integral curve of 
F is up to a time change also a trajectory of the respective Euler flow 
in the upper semi-sphere. For escaping trajectories of F the solution 
of the Euler equation crosses transversally the “equator” {77 = 0} and 
in the lower semi-sphere it follows the same integral curve of F but in 
the reversed direction. Hence it must be a periodic trajectory of the 
Euler flow. 

In the case of a quadratic L — F system we get the following 

Corollary 5.3. If the linear vector field F(£) = Ff has eigenvalues 
both with positive and negative real parts, then for the Euler equation 
(3) the periodic trajectories fill an open and dense subset of E> n C g. 

If all the eigenvalues of F are on the imaginary axis then for the 
Euler equation (3) an open and dense subset of S n C g is filled with 
trajectories which are either periodic or quasi-periodic after an appro¬ 
priate smooth time change. 

Proof. Let us first note that for linear vector field F the condition ( 4 ) 
is satisfied on an open and dense subset of the unit sphere S' C go, 
unless F* = —F. 

If the matrix F has eigenvalues with both positive and negative real 
parts then the set of escaping points in the unit ball B C go is not only 
open but also dense in the unit ball. That is so because the instability 
as f —* Too forces typical solutions to reach the boundary unit sphere 
at some finite time both in the future and in the past. By Theorem 
5.2 they give rise to the periodic solutions of the Euler equations. 
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If the matrix F has only purely imaginary eigenvalues, and there is no 
resonance, then the linear system (5) has only quasi-periodic solutions. 
In general there is an open set of trajectories which stay in the unit 
ball for all times, and an open set of escaping trajectories. The union 
of these open sets is dense in the unit sphere. The only exception 
is the case of the skew-symmetric matrix F* = —F, when there are 
no escaping trajectories. Note that while the solution of (5) is quasi- 
periodic the solution of the Euler equation need not be quasi-periodic 
since a time change is involved. 

In the resonant case we have instability for t — y ±oo, and hence there 
is a dense subset of escaping trajectories, both in the cases of purely 
imaginary nonzero eigenvalues and of the zero eigenvalue. By Theorem 
5.2 we get an open and dense subset of the unit sphere S" C g filled 
with periodic solutions of the Euler equations. □ 

6. Robust integrability of L - F systems 

The notion of integrability of a dynamical system has a long his¬ 
tory and several different versions. In hamiltonian dynamics the main 
concept of integrability is associated with the Liouvillc-Arnold theo¬ 
rem, where families of invariant tori carrying quasi-periodic motions 
appear. The issues involved in the general definition of integrability 
in hamiltonian dynamics were explored by Bogoyavlensky, [Bol] , Bo2 , 
and Fasso [Fj. 

We study a family of reversible, non-hamiltonian systems in which 
quasi-periodic motions occur robustly. We introduce the following 
working definition 

Definition 6.1. We call a dynamical system semi-integrable if an open 
subset U of the phase space is filled with invariant tori carrying the 
quasi-periodic motions. If U is also dense in the phase space then the 
dynamical system will be called integrable. 

It turns out that any periodic solution of the Euler equation con¬ 
structed in Theorem 5.2 gives rise to quasi-periodic solutions of the 
L — F system. This is the contents of the following theorems. 

Theorem 6.1. If a vector field F in go has escaping points in the unit 
ball then the L — F system on SM is semi-integrable. 

If the vector field F has a dense set of escaping points in the unit 
ball then the L — F system on SM is integrable. 

For quadratic L — F system we can formulate an effective criterion 
of integrablity. Moreover the integrablity persists under small pertur¬ 
bations in the space of L — F systems. 

Theorem 6.2. If for a linear vector field F = Ff,f e go the matrix 
F has eigenvalues with both positive and negative real parts then the 
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quadratic L — F system on TM is integrable. Moreover the invariant 
tori are common level sets of real-analytic first integrals. 

Any small perturbation of such a quadratic integrable L — F system 
in the space of all L — F systems must be integrable. 

The last theorem covers the examples of geodesic flows in Butler, [ Bui] , 
and Bolsinov and Taimanov, EH- 

Corollary 6.3. If the operator L : g 0 ~^ flo has eigenvalues with both 
positive and negative real parts then for the left-invariant metric on M 
the geodesic flow on TM is integrable. 

Moreover the (n + 1)-dimensional invariant tori are common level 
sets of real-analytic first integrals in involution. 


It seems that within the class of quadratic integrable systems, for 
an open dense subset of operators F the system is non-degenerate, in 
the sense of having a rich family of frequencies on SM. The reversible 
version of the KAM theory, developed by Moser, [Mj, and Sevryuk, [S], 
would be applicable upon the establishment of the non-degeneracy. The 
allowed perturbations would be J-reversible second order equations, 
but not necessarily left-invariant. 

The calculations required to establish the non-degeneracy are cum¬ 
bersome, and we did not find a satisfactory way to do it. 

We will give a joint proof of Theorems 6.1 and 6.2. 


Proof. We will establish that the periodic solution of the Euler equation 
in g constructed in Theorem 5.2 is covered by invariant tori in TM 
carrying quasi-periodic solutions. To that end, given the T-perioclic 
solution v(t) G Q,v(t) = (£(t),rj(t)) of Theorem 5.2, T = 2 1 0 , let us 
integrate the equations (2) of the L — F system. We get 

u(t) — a + ao(t), ao(t) = / r/(s)ds , 


w(t) — b + b 0 (t), b 0 (t) — f e u ^ L f(s)ds. 

Jo 

Since the function r/ is odd and periodic we obtain that ao(t) is also a 
periodic even function and assumes values between — w 0 and u o, where 
u 0 = /q° r]{s)ds. Note that although b 0 (t) is not periodic, its time 
derivative is T-periodic. 

Fixing the value a = w(0) we obtain an invariant torus. Indeed con¬ 
sider the flow x I ,f defined by the L — F system (3), which is a group 
extension flow by Proposition 3.1. The mapping T t takes M x {u(0)} 
into itself and it acts there as a right translation. Since u{t) is periodic 
the translation is by an element from Go, and hence it preserves all 
the toral leaves. For a G R and v G g denote by £(a,v) the toral 
leaf £(a,v) = {(uyu; f,7]) G SM\w mod To,w = = r;}. Since 

T 7 takes £(a,r>(0)) into itself, and it acts there by the translation 
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by Co = Jq* 0 e u ^ L ^(s)ds, then the union of the n -dimensional tori 
\h*/C(a,u(0)), t € M, is an (n + l)-dimensional torus carrying the sus¬ 
pension flow of the toral translation, i.e., a quasi-periodic flow. 

In the case of a linear vector held F(£) = F£, for some constant 
matrix F, the real analytic hrst integrals are given by the vector valued 
function 

$(w,u;£,ri) = e~ uF £. 

These hrst integrals are functionally independent, hence their common 
level sets must be in general unions of submanifolds of dimension n +1. 
Clearly the toral leaves £(a,v) belong to the level sets, and so do the 
n + 1-dimensional tori constructed above. Hence in general the level 
sets must be unions of the (n + l)-dimensional invariant tori. □ 

The geodesic how on TM is both an L — F system and hamiltonian 
but it is an exception. The L — F systems are rarely hamiltonian. 

To establish the Corollary 6.3 we consider more generally hamilton¬ 
ian systems on T*M with Go symmetry. If the hamiltonian function 
H : T*M —y R is invariant under the action of the torus Go/T 0 then 
we get immediately n hrst integrals in involution. Indeed, this action 
in the canonical variables ( w,u;p w ,p u ), associated with the variables 
(w, u ) in M, amounts to translations in w with the other variables hxed. 
It follows that the hamiltonian FI does not depend on the variables w 
and so the n momenta p w are hrst integrals of our hamiltonian system. 
Together with H itself we have n + 1 hrst integrals in involution, and if 
H is functionally independent of p w then the Arnold-Liouville Theorem 
is applicable. In particular we obtain invariant tori on compact level 
sets of H. 

Let us identify the cotangent bundle T*M with the tangent bun¬ 
dle TM using the left invariant metric on M. This identification is 
the Legendre transform associated with the geodesic how. Since the 
lagrangian L for the geodesic how is equal to 

JC(w, u ; w, u) — ^ (( e~ uL w , e~ uL w) + u 2 ) , 
we get immediately that 

p w = vprL = e~ uL * e~ uL w = e~ uL * f. 
ow 

For the geodesic how these are the same hrst integrals as those in 
Theorem 6.2. 

Let us finally consider more special hamiltonians, which are invariant 
under the action of the full group G. They can be described by func¬ 
tions H : g —> R, H — H(£, T)) and H = H(e uL *p w ,p u ). It is instructive 
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to compare the respective hamiltonian equations in the (w, u ; £, rf) vari¬ 
ables in TM with the equations of an L — F system. We have 

= ( h_ = fdH 

dt di] ’ dt df 
dw uL dH du dH 
dt df ’ dt dr) 

It transpires that also in the hamiltonian case if the function H has 
compact level sets then the Euler equation in g has an open and dense 
set of periodic trajectories. 


7. Compact configuration spaces 

In this section we consider the additional properties of an L — F 
system on the compact phase space, namely the unit sphere bundle 
SN. where T is a cocompact lattice in G and N = T\G, as discussed 
in Section 1. 

Let us note first that in this case the group must be unimodular, 
and the operator L has zero trace. Hence either L has eigenvalues with 
both positive and negative real parts, or all of its eigenvalues are on 
the imaginary axis. 

The integrability and semi-integrability are passed from the non¬ 
compact phase space SM to the compact phase space SN without 
further assumptions. What is new is the appearance of hyperbolic 
behavior made possible by the recurrence in the compact phase space. 
For any quadratic L — F system the submanifold A C SM given by the 
equations £ = 0, r\ — 1 is diffeomorphic to M and it carries trajectories 
escaping to infinity both in the future and in the past. The projection 
of A to the compact phase space SN is diffeomorphic to N and the 
L — F flow is the suspension of the toral automorphism e L . If the toral 
automorphism has no eigenvalues on the unit circle then it is an Anosov 
diffeomorphism, and its suspension is an Anosov flow. If it has only 
some eigenvalues outside the unit circle then we get partially hyperbolic 
flows as suspensions, im Such flows in the compact phase space 
SN while integrable have positive topological entropy. The existence 
of such flows was the discovery of Bolsinov and Taimanov, [eh. 

More generally let us assume that the automorphism e L has eigen¬ 
values outside of the unit circle. Then for an arbitrary left-invariant 
second order equations on SN (not necessarily an L — F system) we 
have the following 

Proposition 7.1. If the variable rf has positive lower (upper) time 
average on some trajectory then this trajectory has negative and positive 
lower (upper) Lyapunov exponents. The dimension of the respective 
stable and unstable subspaces is equal to the dimension of the stable 
and unstable subspaces of the automorphism e L . 
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Proof. The flow has the structure of a group extension flow over 
the Euler flow as described in Proposition 3.1. We fix a trajectory 
v(t) = (£o(£), Vo(t)) °f the Euler flow. There is a unique parametrized 
curve (w 0 (t),u 0 (t)) G G,U7 0 (0) = 0,u 0 (0) = 0, such that for every 
(w, u) G G 

T^-uyu; u(0)) = (w + e uL w 0 (t), u + u 0 (t); v(t)) . 

Hence the flow takes the n- dimensional toral leaf £(rqu(0)) onto the 
toral leaf L{u + u 0 ,v(t)). We have the coordinates w for both toral 
leaves. In these coordinates the restriction of the flow is defined by 
w i—>■ w + e uL Wo(t). Differentiating this mapping with respect to w we 
obtain that Dfid restricted to the n- dimensional tangent subspace of 
the toral leaf is given by the identity operator in the coordinates w. 
However we measure the tangent vectors using the left-invariant metric 
and that means 

\\dw\\ 2 = (e~( u+uo)L dw) 2 . 

Hence we get the exponential growth (decay) as long as u 0 has linear 
growth (decay). Since = ho, the assumption of the positive time 
average of r) 0 leads to the exponential growth (decay) of ||chn|| 2 in the 
unstable (stable) subspace of the operator L. □ 

We have not chosen any invariant measure so the Lyapunov expo¬ 
nents need to be understood as the upper (or lower) limits. It follows 
from the proof that assuming negative lower (upper) time average we 
will arrive at the same conclusion, with the reversal of the number of 
positive and negative Lyapunov exponents. 

The phenomenon described in Proposition 7.1 appears explicitly in 
the work of Butler and Paternain on special magnetic flows on SOL, 
[BmP] , 

Let us consider first integrals of a quadratic L — F flow HP on the 
compact phase space SN. The real-analytic first integrals of Theorem 
6.1 can be considered as real analytic multi-valued first integrals. We 
will establish that for a residual subset of integrable quadratic L — F 
systems on SN there are no single valued real-analytic first integrals 
constant on the toral leaves. If an integrable system is non-degenerate 
(i.e., it has a rich family of frequencies for the quasi-periodic motions 
on the tori) then by necessity any continuous first integral must be 
constant on the toral leaves. Since we have not established the non¬ 
degeneracy we need to include the constancy into the assumptions. 

Let us recall that by Theorem 6.1 quadratic L — F flows on SN 
for which the matrix F has eigenvalues with positive and negative real 
parts are integrable, and they form an open subset in the space of all 
quadratic L — F flows. The complement of this subset has an interior 
filled with flows asymptotic to the suspension of the toral automor¬ 
phism e L . 
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Proposition 7.2. For a residual subset of integrable quadratic L — 
F flows on SN there are no single valued real-analytic first integrals 
constant on the n-dimensional toral leaves, i.e., functions of (u]f,r]). 


Proof. We consider only diagonalizable operators F, with all different 
eigenvalues (Ai, A 2 ,..., A n ), at least some of them with positive and 
some of them with negative real parts. This gives us an open and 
dense subset of integrable quadratic L — F systems. If there is a real- 
analytic first integral R for such a system on SN then it can be lifted to 
the first integral of the system on SM, which we denote also by R. Now 
R can be considered as a function of (u; £, rf). We restrict our attention 
to the unit ball B in g 0 (the “£-space”), and we put 77 = a/1 — fl 2 . This 
gives us a real-analytic function 

R{u\ £), ( 77 ; £) G R x B, which is also a first integral of the system. 

We further consider the time change ^ = 77 in our domain 1x5 
The equations of the system in the variables ( 77 ; £) are 


( 6 ) 


— = 1 . ^=F£ 

ds ’ ds 


The function R is a first integral of ( 6 ) so that its time derivative is 
identically zero 


DR 

du 


= 0. 


It follows that for the real analytic function P{u\Cf) = R(u;e uF £ ) de¬ 
fined in the vicinity of the origin we have 


dP 

du 


dR dR 
du df 


Fe uF C = 0. 


Hence P is the function of f alone defined in the vicinity of the origin 
in R n such that R(u;f) = P(e~ uF f). Now we use the fact that the 
function R(u\£) = P(e~ uF f) must be actually a periodic function of u 
with the period 1, since R is lifted from SN. Hence 


(7) P(e~ F f) = P(0, for all £ G R n 


It goes back to Poincare that such a relation is impossible if the eigen¬ 
values of F are non-resonant. More precisely, if there is no integer 
vector (fci, k 2 ,..., k n ), with non-negative entries, such that 

k\ Ai + - • -+k n X n = 0 then the relation (7) can hold only for a constant 
function. 

The family of L — F systems with non-resonant eigenvalues form a 
residual subset of the open set of integrable quadratic L — F systems 
introduced at the beginning of the proof. □ 
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8. Invariant measures 

In this section we discuss smooth invariant measures for the L — F 
systems, for the non-compact SM and compact SN phase spaces. First 
of all we consider the Lebesgue measure dw du in G, which is the right- 
invariant Haar measure. It is also invariant under the left translations 
by elements from To. So it projects into a a-finite measure p on M = 
T 0 \G. The product of p by the standard Lebesgue measure in § n C g 
(or g) will be referred to as the Lebesgue measure in SM — M x 
(or TM ) and denoted by v. Once the Lebesgue measure v is chosen 
we call a function p > 0 an invariant density of a dynamical system if 
the measure pv is preserved by the dynamical system. 

Proposition 8.1. If the vector field F has constant divergence in the 
unit ball of g 0 then p = exp(—u div F) is an invariant density for the 
L — F system on SM. 

Proof. Since the L — F system defined in TM — M x g preserves all 
the sphere bundles it suffices to check the invariance of the measure in 
TM. It can be easily seen that the divergence of the L — F system in 
TM is equal to p div F. Since 

—p = —77 div F, 
dt 

the claim is proven. □ 

In the case of a compact phase space SN the group G must be 
unimodular, i.e., tr L = 0. For unimodular Lie groups the left Haar 
measure is equal to the right Haar measure, and the Lebesgue measures 
/x and v project as finite measures to N and SN, respectively. We will 
denote the resulting measures again as p and u. 

It follows from the above Proposition that if div F = 0 then the 
L — F system preserves the Lebesgue measure v in SN. Moreover the 
Euler flow of the L — F system preserves the Lebesgue measure on 

S n C g. 

Conversely, if an L — F system has a (real-analytic) invariant density 
p in SN then the Euler flow in C g has a (real-analytic) invariant 
density. Indeed, we obtain the density of the projected measure by 
integrating out the w, u coordinates. It preserves the smoothness, and 
the real-analyticity of the density. 

For the quadratic L — F systems the Euler flow cannot have a smooth 
invariant density which is positive at the fixed point f = 0 , unless 
div F = 0. However there may be smooth, or even real-analytic, in¬ 
variant densities which vanish at the fixed point, even if div F^O. 

Proposition 8.2. For a residual subset in the space of quadratic L — F 
systems in SN there are no real-analytic invariant densities. 

For the open family of integrable quadratic L — F systems in SN 
with the matrix F having all different real eigenvalues, there are C°° 
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invariant densities. Moreover there is a dense subset in this family with 
real-analytic invariant densities. 


Proof. It follows from the above considerations that if an L — F system 
in SN has a real analytic invariant density then the field F has a real 
analytic invariant density p = p(£) in the unit ball in g 0 . Hence we 
have 

div ( pF) = 0. 

For a diagonalizable matrix F we can analyze the last equation in 
the coordinates in which F is diagonal. This may require a linear 
complex change of coordinates, but that is not a problem in calculations 
involving Taylor expansions. Assuming that F is already diagonal with 
eigenvalues Ai, A 2 ,..., A n we get 


( 8 ) 


i= 1 


dp 

dii 


—trFp 


This relation can hold for a function real-analytic in the neighborhood 
of the origin £ = 0 if and only if there are monomials which satisfy it. 
For p(£) = niLr £ ri we get substituting into (8) 


(9) A j{rj + 1) — 0. 

i —1 


The set of matrices F for which (9) does not hold for any vector 
(rq, r 2 ,..., r n ) with natural entries is a residual subset in the space 
of all matrices. The first part of the Theorem is proven. 

To prove the second part consider the set of matrices with all different 
real eigenvalues, at least one positive and one negative. It is an open set 
of matrices. Its subset where the eigenvalues satisfy (9) for some vector 
(ri, rq,..., r n ) with natural even entries is dense. For such systems the 
function p(£) = Y\ff =l f ri > 0 is an invariant density of the Euler flow, 
and also the full L — F system in SN. 

To get a C°° invariant density we choose a vector (r 1; r 2 ,..., r n ) 
with positive entries satisfying (9) and consider the function p(£) = 
mu ir > 0. The function satisfies (8) in the open dense subset 
where it is positive, but it is not in general a smooth function. We 
consider the first integral of the Euler flow /(£) = n”=i |£| ri+1 an d we 
get a C°° invariant density p(£) = exp(— /(£) _1 )p(0- CD 


The C°° invariant densities in the proof are similar to the C°° first 
integrals of Butler, Bui . 


9. Coexistence of integrablility and positive metric 

ENTROPY 

The results of Sections 6 and 7 may be used to find real-analytic vol¬ 
ume preserving L—F flows with both quasi-periodic motions in an open 
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subset and positive metric entropy. Such examples were constructed 
by Butler in the C°° class, [Bu2] . [Bu3 ]. The recent survey of Chen, Hu 
and Pesin, (C-H-Pj describes other kinds of coexistence phenomenae. 

Let F have purely imaginary all different eigenvalues, but F* —F. 
For such a linear vector held the open unit ball is not invariant, hence 
there is an open subset of escaping points. More precisely the whole 
open unit ball is the closure of the union of two open subsets V e and 14. 
The open subset V e contains only escaping points, and the open subset 
14 contains points whose trajectories have compact closures in the unit 
ball. Actually these compact closures are invariant tori of the linear 
system defined by the matrix F. Note also that the linear vector held 
has zero divergence, so the L — F how preserves the Lebesgue measure. 

These considerations apply already for n = 2. If A has purely imag¬ 
inary eigenvalues, but F* —F, then the open unit disk contains an 
open ellipse filled with elliptical integral curves, The rest of the open 
unit disk is filled with escaping trajectories. 

The open set of escaping trajectories V e gives rise by Theorem 6.1 to 
semi-integrability of the L — F how. If the matrix L has eigenvalues 
with positive (and negative) real parts then the open set of bounded 
trajectories 14 leads by Proposition 7.1 to positive Lyapunov exponents. 
By Pesin formula, ra, we obtain also positive metric entropy with 
respect to the invariant Lebesgue measure in SN. We get 

Theorem 9.1. There are volume-preserving quadratic L — F systems 
on SN which are semi-integrable and have positive metric entropy. 

Let us note that these properties will hold also for small divergence 
free perturbations of F as long as we guarantee that the perturbed 
vector held in g still has a positive Lebesgue measure of quasi-periodic 
motions near the origin £ = 0. This can be achieved for instance 
when the linear vector held F and the perturbation are hamiltionian 
and satisfy the non-degeneracy conditions of the KAM theory, (A], 
Appendix 8. 

Further let us consider the modified linear vector held Fi = F — el. 
We obtain for e > 0 the asymptotic stability of the origin £ = 0. At 
the same time for small e the restriction of the vector held Fj to the 
unit ball still has an open set of escaping trajectories. 

Theorem 9.2. There is an open set of quadratic L — F systems on 
SN which are semi-integrable and the interior of the complement to 
the set of invariant tori is filled with orbits which are asymptotic as 
t —> +00 (—00) to the attractor (repellor) carrying the suspension of 
the toral automorphism e L (e~ L ). 

We conjecture that the following is true. 

Conjecture. There are real-analytic divergence free vector fields in 
the closed unit ball with an open and dense subset filled with escaping 
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trajectories, and with the complement of positive Lebesgue measure, 
filled with trajectories defined for all times and which are contained in 
the interior of the unit ball. 

In the Arnold diffusion scenario a generic hamiltonian perturbation 
of an integrable hamiltonian system has a nowhere dense subset of 
invariant tori of positive Lebesgue measure. In the complement almost 
all orbits move away unboundedly both in the future and in the past. 
However this scenario was not so far rigorously established in any real- 
analytic examples. 

If this conjecture holds then the respective real-analytic L — F system 
would be integrable and of positive metric entropy. 

10. Gaussian thermostats 

A Gaussian thermostat is defined by a Riemannian metric on a man¬ 
ifold M and a vector field E. The trajectories of the Gaussian thermo¬ 
stat satisfy the following ordinary differential equations in the tangent 
bundle TM, (EWj. 

(10) x = v, V> — E — 

dt (v , v) 

where x = x(t) G M is a parametrized curve in M. 

By the force of these equations the “kinetic energy” v 2 is constant. 
Fixing the value of this constant k = v 2 , and introducing the auxiliary 
vector field F\ = j^E we can rewrite the equations (10) as 

(11) ^ x = v, V v v = v 2 F k - (F k , v)v. 

It was observed in [P-W] that the equations (11) describe geodesics of 
a special metric connection V defined by the field F k 

S7 X Y = S7 X Y - (X, Y)F k + (Y, F k )X, 

where A", Y are arbitrary smooth vector fields on M. This connection 
is not symmetric, it is the unique metric connection with the torsion 
T(X, Y) = (Y, F k )X - (X, F k )Y, [P=W]. 

Note that the equations of a Gaussian thermostat define significantly 
different flows for different values of k — v 2 , while the equations (11) 
scale so that for different values of v 2 we get the same flow up to a 
change of time by a constant factor. To get complete understanding of 
the dynamics of (10) we need to consider the dynamics of the geodesic 
flow of the metric connection (11) for the whole family of vector fields 
F k = \E, k > 0. 

If the manifold M is the locally homogeneous space of Section 1 
and the vector field E is left invariant, then we obtain a family of 
left-invariant para-metric connections. To get a weakly 3C-invariant 
connection we have to have KE = E, i.e., E must be orthogonal to 
0o- The equations of geodesics for these connections give us a family of 
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quadratic L — F systems. Assuming that E 2 = 1, we get that in this 
case the respective vector held in g 0 is equal to 

F(0 = 

We obtain the following corollary of Theorem 6.2 and the considerations 
leading to Theorem 9.2. 

Corollary 10.1. Let r min and r max be the smallest and the largest real 
parts of the eigenvalues of L. 

If r min < r max then for k = v 2 such that r min < \ < r max the 
Gaussian thermostat (10) in M is integrable. 

For k = v 2 such that I > r max the trajectories of the Gaussian 
thermostat (10) in the compact space N are asymptotic to a subsystem, 
which is the suspension of the automorphism e L . 

In particular if r min < 0 < r max then for 0 < k < —— the Gaussian 
thermostat (10), for the held E orthogonal to go, has a global attractor 
carrying the suspension of e L (which can be an Anosov how). However 
for k > — — it is integrable. This is a somewhat paradoxical behav- 
ior: strong mixing properties at low kinetic energy, and quasi-periodic 
motions for high kinetic energy. Admittedly the strong mixing proper¬ 
ties occur at the asymptotic submanifold of half the dimension of the 
phase space. And this subsystem is present for all values of the kinetic 
energy, but it is not an attractor or a repellor in the integrable case. 


11. A FAMILY OF RlEMANNIAN METRICS WITH INTEGRABLE 

GEODESIC FLOWS 


The Gaussian thermostats can be associated with Weyl connections, 
|Wlj . |W2j . The left invariant vector held E of Section 10 is a gradient 
of a function on the group G, which factors to M — T 0 \G. Indeed we 
have E = Jjj = Vu. The Weyl connection defined by gradient vector 
helds are Levi-Civita connections of modified metrics, namely for the 
vector held Fk = jE we need to multiply the initial left invariant metric 
by the function eir. 

Let for simplicity L be diagonal with eigenvalues Ai, A 2 ,..., A n . The 
left invariant metric on G is given by 

ds 2 = e~ 2uXl dw\ + e~ 2uX2 dw\ H-b e~ 2uXn dw 2 n + du 2 . 


Introducing the new variables 


x 0 = e>=,Xi = T Wi,i = 1, 
k 


we obtain 


e 2 k ds 2 = k 2 (^x 2 Q 1 dx\ + x(f 2 dx 2 2 + • • • + x(f' n dx 2 n + dxl) , 

where t* = 1 — k\i,i = 1 ,n. It was calculated in jT-Wj that the 
Weyl sectional curvatures are non-positive in this case if and only if 
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Tj > 1 or r, = 0 for i = 1 ,n. The sign of sectional curvatures of 
the last metric must be the same as for the Weyl connection ( [P-Wj h 
Considering the geodesic flow of the Riemannian metric on M = 
T n x M + = {(cc; a; 0 )|cc mod r 0 ,x 0 > 0}, where T 0 is a lattice in R n , 
we conclude that it is integrable if and only if there are exponents r* 
of opposite signs. It follows from the considerations of Section 10, but 
this phenomenon can be greatly generalized. Note that the metric is 
not complete since the variable xq is assumed positive. Nevertheless 
in the integrable case geodesics from an open and dense subset can 
be extended indefinitely since they do not leave a compact subset of 
M. If the conditions of integrability are not satisfied then almost all 
geodesics leave every compact subset of M. 

For an open interval (a, b ) C R, finite or infinite, let us consider a 
Riemannian metric on 

M — T n x (a, b ) = {(x; xo)|x mod r 0 , a < x o < b}, 

ds — oqcLtq T cx^dx^ T ■ ■ ■ T ct n dx n T dx g, 

where all the functions cq are positive functions of the variable xq alone, 
defined on the interval (a, b). 

Proposition 11.1. If lim Xo _^. a aflxo) = 0 and linn^^f, cq(xo) = 0 for 
some i and l then the geodesic flow is integrable. 

Proof. Let us pass to the hamiltonian formulation. We get the mo¬ 
menta 

Po = X 0 , Pj = a)xj,j = 1 ,..., n, 
and the hamiltonian 

n =\ +it a 7 2 ^ • 

We have the obvious n first integrals pi,...,p n . Together with H we 
get n+1 first integrals in involution. Under the assumption that Pi 0 
and Pi 0 the common level set of the first integrals must be compact 
in the variable x 0 , hence altogether compact. Indeed, we have 

H > oif 2 pl + otf 2 p 2 —> +oo as xq —y a, b. 

By the Arnold-Liouville Theorem all regular compact level sets are tori 
carrying the quasi-periodic motions. □ 

This last Proposition takes us away from homogeneous systems, now 
we need only the T" symmetry. 
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